log P_n: Calculate marginal probability of observed genotypes under 'null' model
Description
Marginal probability calculated exactly by integration.
Usage
"log P_n"(y, q_shape = c(2, 100))
Arguments
y
Logical vector of case (TRUE) control (FALSE) status.
q_shape
Beta shape hyper-priors for prior on rate of affection (i.e. being a case) amongst individuals with non-pathogenic variant combinations (i.e. they have less than min_ac variants.